AN OPTIMIZATION PROBLEM WITH FREE BOUNDARY 
GOVERNED BY A DEGENERATE QUASILINEAR OPERATOR. 



KRERLEY OLIVEIRA AND EDUARDO V. TEIXEIRA 



Abstract. In this paper we study the existence, regularity and geometric 
properties of an optimal configuration to a free boundary optimization problem 
governed by the p-Laplacian. 



1. Introduction 



Let D be a bounded domain in R n and tp a positive function defined on it (the 
temperature distribution of the body D). A classical minimization problem in 
heat conduction asks for the best way of insulating the body D, with a prescribed 
amount of insulating material in a stationary situation. This situation also models 
problems in electrostatic, potential flow in fluid mechanics among others. 

The mathematical description of this problem is as follows: fixed a number, 
7 > (quantity of insulation material), for each domain f2 surrounding D, such 
that |f2\£)| = 7, we consider the potential u associated to the configuration f2, i.e., 
the harmonic function in il\D, taking boundary data equal to tp on dD and on 
d£l. The flow of heat (quantity to be minimized) corresponding to the configuration 
Q, is given by a nonlocal monotone operator given by, 



where [i is the inward normal vector defined on dD. The function T: dD xM^M 
is assumed to be convex and increasing, on and continuous on x. Important 
examples are, T{t) — t (classical heat conduction problem), T(t) = t p (optimal 
configurations in electrostatics), T(x,t) = max{i, C(x)} (problems is the material 
sciences) . 

In |Teix05j . the second author studied the problem of minimizing J(f2) among 
all configurations such that, say, |f2 \ D\ = 1, where \A\ is the volume of the set 
A. This optimization problem with linear heat flux, i.e., T(x,t) = t, was studied 
in AAC86] and |AUS88| . Qualitative geometric properties of the free boundary, 
namely symmetry, uniqueness and full regularity of the free boundary, was explored 
in |Teix-lj . 

In this present paper, we turn our attention to this problem when we allow the 
temperature itself deform the medium. We assume the influence of the temperature 
distribution on the medium is proportional to the magnitude of its gradient. These 
considerations lead us to study this optimization problem when temperature dis- 
tribution is governed by the p-Laplacian. In other words, the variational problem 
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we are interested in is, for 1 < p < oo 



(1.1) Minimize 



J(u) := T(x,u li )dS u: D c -> E, u = ip on 3D, 

JdD 

A p u = in {u > 0} and \{u > 0}| = 1 



Here, A p u := div(| Vw| p ~ 2 Vu) . Several new difficulties appear when dealing 
with the nonlinear operator A p . One of its main difficult lies in the fact that the 
p-Laplacian is not uniformly elliptic. 

In analogy with the linear case for the Laplacian operator, in this paper we 
shall restrict ourselves to the heat flux given by T(x,t) = t p_1 . In our physical 
considerations, we will assume the body to be insulated has much smaller volume 
than the quantity of insulation material. This leads us to consider a constant 
temperature distribution, say, ip = 1. With these assumptions, problem (|1 . If) can 
be reformulated in terms of the following equivalent version of it: 

u: D c — » M, u = ip on 3D, 
A p u = in {u > 0} and \{u > 0}| = 1 



...... , J(u) := / \Vu\ p dx 

(1.2) Minimize { K ' J D c 1 



In problem (|1.2|l , we may assume tp to be a mere continuous and positive function. 
When ip = 1 in problem l|1.2|l . such a problem is equivalent to problem 1)1. l|l with 
T(x, t) = fP^ 1 . In this paper we will only deal with problem (|1.2I) . We hope to turn 
our attention to problem 11.1(1 in its full generality in future research. 

From the mathematical point of view, our approach is motivated by recent ad- 
vances on the free boundary regularity theory for minimum problems with a vari- 
able domain of integration involving degenerate quasilincar operators. Namely, D. 
Danielli and A. Petrosyan, in | DP05| . have recently extended the celebrated work of 
H. Alt and L. Caffarelli AQ2], for the p-Laplacian operator. Our regularity results 
will relay on suitable modifications of the arguments in |DP05| . Furthermore, we 
shall establish a free boundary condition that will relate our optimization problem 
with Bernoulli type problems, similar to the ones studied in [HSfT)] . jDPS03 and 

[DM!. 



2. Mathematical fundaments of the physical problem 

In this section we shall introduce the main mathematical tools we shall use 
throughout the whole paper. Throughout the article, 1 < p < oo and A p u stands 
for the p-Laplacian operator 

A p u := div(|V|^ 2 Vu). 

Let U be a domain in R N , Let us recall that for any £ G W 1 '?^), A p £ G 
[W£' P (U)]* and 

{A p £,v)= [ \V£,\ p - 2 V£,-\7vdx. 
Ju 

Problem (|1.2|) presents several difficulties from the mathematical point of view. 
Our strategy will be to study a penalized version of it, which is sort of a weak 
formulation of the problem. The idea is to grapple with the difficulty of volume 
constraint, which is very unstable under limits and makes perturbation arguments 
quite hard. 
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From now on, we denote by V the following set 

V := ju G W 1 ' P {D C ) u = cp on <9L>} . 

The penalized problem is stated as follows: Let e > be fixed. We consider the 
function 

, _ / 1 + ~(t- 1) if t > 1 

[ 1 + e(t — 1) otherwise. 
We then define the penalized functional as 

(2.1) J s (u) := J \Vu\Pdx + f s (\{u > 0}|). 

For the moment, we shall be interested in the following minimization problem 

(2.2) minJ e (u). 

For latter use, given a 7-bilipschitz function / we consider the similar penalized 
problem of l|2.2|) given by the functional: 

(2.3) Jf.= [ \Vu\ p dx + f(\{u > 0}| 

J DC V 

3. Properties of solutions of Problem ()2.2|) 

In this section we shall derive existence of a minimizer for the penalized problem 
as well as some important nondegeneracy conditions, such as optimal regularity 
and linear growth away from the free boundary. The proof of Theorem 13.11 will 
be developed throughout this section. At the end, we shall be able to state a 
representation theorem that will be crucial to study further regularity properties of 
the free boundary. 

Theorem 3.1. For each e > fixed, there exists a minimizer u t £ V for the 
functional J E . Furthermore 

(1) u e > 0. 

(2) A p u e is a nonnegative Radon measure supported on d{u e > 0}. In partic- 
ular, 

A p u E = in {x G D c u £ (x) > 0}. 

(3) Any minimizer u e of problem is Lipschitz continuous and for any 
compact K, C D there exists a constant K — K(K.,p,n, e, D) such that 

||Ue||itp(/C) < K. 

(4) The function u e growths linearly way from the free boundary, i.e., for any 
compact K, C D c , there exist positive constants c, C, depending on dimen- 
sion, IC, p, D and e, such that 

cdist{x,d{u E > 0}) < u e (x) < Cdist(x,d{u e > 0}), Vx € IC. 

(5) The free boundary is uniformly dense, i.e., for any compact JC C D c fixed, 
there exist a constant c = c(IC,p, n, e, D) with < c < ^> such that 

\B r n{u>o}\ 

\B r \ C ' 

for any ball B r = B r (x) centered at some point x G d{u e > 0} PI K,. 
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Proof. The existence of a minimizer u £ for problem (|2.2|1 follows easily from the 
fact that for any minimizing sequence u k , we may assume 

• Vu k -± Vm e in LP(D C ) 

• u k — > u £ almost everywhere in D , 

for some u £ G V. Thus, 

\\7u £ \ p dx < liminf / \X7u k \ p dx 

D C fc-K» J D c 

and 

> 0}| < liminf \{u k e > 0}|. 

k — >oo 

Since / e is a continuous and increasing function, we obtain 

Js(u £ ) < liminf J{u k e ). 

k — >oc 

Clearly, u £ > 0, otherwise, J £ { (u £ ) + ) < J e (u £ ). 

Observe that if u £ is a minimizer, then J{u £ ) < J(u £ — erj), for every e > and 
nonnegative rj g Cfi°(D c ). Since {u E > erf} C {ue > 0} and / is increasing, we 
have that f({u £ > 577}) < /({u e > 0}) and consequently, we have that A p u £ is a 
nonnegative Radon measure supported on d{u £ > 0}. 

Now, we explain the main track and the necessary changes in order to obtain 
items (3), (4) and (5) . We will follow the lines of [DPQ5j establishing a sequence of 
lemmas that are analogous of those in this refereed work. Instead of proving these 
lemmas with all the details, we shall restrict ourselves to enunciate them and sketch 
their proofs by pointing out the necessary modifications and sending to DP05 for 
further details. 

Up to the end of this section, we fix some j-bilipschitz function f . Consider u 
a minimizer of the problem 12.31 in some ball B. Denote by v the solution of the 
Dirichlct problem 



(3.1) 



A p v = in B 



v = u on dB. 
Following the beginning of Section 3 in DP05] we obtain: 
Lemma 3.2. There exists a constant C = C(n,p, 7) such that 



|V(w -v)\ p < C\{u > 0}n B\ for p> 2 and 

( 3 - 2 ) r JB r 

/ |V(u- v)\ p < C\{u > 0}nS|5( / \Vu\p) 1 1 forl<p<2. 

J B J B 

Moreover, the constant C goes to zero when 7 goes to zero. 

The next Lemma is the analogous of Lemma 3.1 in DPO^j, with a similar proof. 

Lemma 3.3. Let f be a given 7- Lipschitz function and u be a bounded minimizer 
of Problem \2.!ft in B\. Then, u € C' a in Br for some a = a(n,p) G (0, 1) and 

\\u\\ C a(B 7 ) < C{n,p, ||u|| L oo (Bl) , 7 ). 
Following Lemma 3.2 in |DP05| . we obtain 
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Lemma 3.4. Let u be a bounded local minimizer of Problem ] 2., -A in B\ with u(0) 
0. Then, there exists a constant C — C(n,p, 7) such that 



3 



< c. 



Proof. Assume, by contradiction, there exists a sequence Uk of bounded local min- 
imizcrs of Problem 12.31 in B\ with u k (0) = and max^j u k > k. In this case, we 

4 

define: 

d k (x) = d(x,d{u k > 0}) 

and 

O k = ^xGB 1 ;d k (x)<^^ 

Observe that, since u k (Q) = 0, d k (x) < \\x\\ for every x G B\. On the other hand, 
(1 — || x ||)/3 > 1/4 for every x G B1/4. From this inequality we conclude that 
S1/4 C O k . Now, define 

3 3 
m k := max(l — ||a;||)itfcfx) > - maxu k (x) > —k. 
o k 4 s 1/4 4 

Consider any maximum point x k G O k of the function (1 — ||x||)'Ufc(a;) and observe 
that 

(3.3) u k (x k ) = - _ > -k. 



Denote by y k any point in d{u k > 0} such that d k {x k ) = \\y k — x k \\ and define 

3 



<5fe ;= \\ x k — Vk\\- Since x k G we have that 5 k < — jp^ . Then, for every 



z e B26 k {yk)- 

\\z\\ < \\y k \\ + 24 < K|| + 34 < KH + (1 - KID < 1. 

From this inequality, we have that B2& k (y k ) C B\. Now, we claim that B$ k /2(y k ) C 
O k . In fact, if z G Bs k /2(yk) we have that: 

4 < < (iHW!) 

w - 2 - 3 - 3 

where the second inequality is a consequence of ||j/fe|| + 4 < !■ Moreover, for 
z G B s/2 (y k ) 

(1 - INI) > (1 - KID - Hi* - *|| > (1 - \\x k \\) -ls k > {{ " ' ! 



2 2 

As a consequence of this inequality, 

(3.4) _max u k < 2u k (x k ). 

By the definition of S k , we have that Bs k (x k ) C {u k > 0}. Recall from item (2) of 
Theorem 13. II that A p u k — in B$ k (x k ). By the Harnack Inequality for p-harmonic 
functions, we may conclude that there exists a constant c = c(n,p) > such that 

(3.5) _ min u k >cu k (x k ). 
In particular, 

(3.6) _max u k > cu k {x k ). 
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Consider the following scaling of Uk 

(3.7) Wk(x) = — f — , for x G Bi 

u k (x k ) 

Observe that since a is a local minimizer of Problem 12.31 Wk is a minimizer of 
the analogous problem replacing f £ by //. = f e /uk{%k)- I n other words, Wk is a 
local minimizer of 

j fc H= / ||v w p + / fc (|{ w >o}|). 

•'■B3/4 

Now, we denote by «fc the solution of A p Vk = in B 3 u and Vk — Wk G Wj' p (B 3 / 4 ). 
Initially we observe that ||/fe||z,i P converges uniformly to zero, as k — > oo. Since u>& 
is a minimizer of , by Lemma 13.21 we guarantee the existence of a sequence of 
positive numbers Cfc, that converges to zero as k — > oo, such that 

(3.8) / \\V(w k -v k W<C k . 

By estimates l|3.4[) and <|3.6[1 we have that 

msptwk < 2, maxwfc > c and iufc(0) = 0. 

Bi b 1j/2 

Now notice that Wk is uniformly bounded, thus from Lemma 13.31 we conclude w k 
and Vk are uniformly C™ in B 5 / s . By compactness, we may assume (passing to a 
subsequence, if necessary) Wk — > wo and ^fc - * uniformly on B 5 / 8 . Equation l|3.8(l 
implies that wo = vq + K in B 5 / 8 . Since Ufe — ► wo, we have that A p vq = and this 
implies that A p wq = 0. By the strong minimum principle, we have that wq — in 
S5/8, because wo > and wo(0) = 0. It is a contradiction with maxw > c > 0. 

Bl/2 

This finishes the proof. □ 

Now, we show how one can derive Lipschitz continuity of a solution u of Problem 
12.31 using Lemma T3. 31 and Lemma f3. 41 

Proof of Lipschitz continuity of u: First, note that it is sufficient to prove 
that u is Lipschitz continuous at every point y on the boundary of {u > 0}, since 
A p u = on the open set {u > 0}, and therefore u is C lo '" on {u > 0}. Take 
y € d{u > 0} and consider u r : B\ — * K defined by: 

u r (x) := -u(y + rx). 
r 

Since u is a minimizer of Problem l2.3l in the ball B r (y), we obtain it r is a local min- 
imizer of the same problem in the ball B±, with u r (0) = u(y) — 0. By Lemma 13.41 
we have that ||urlU° o (B 1 / 4 ) < where C is a constant depending only on p, 7 and 
n. Thus, by Lemma 13.31 we conclude u r G C a (B 1 / 4 ). Furthermore, for every r > 

||w||g°(b 1/4 ) < C = C{n,p,i). 

This implies that ||Vu||oo < C and thus u is locally C-Lipschitz (see Theorem 4.2.3 
in |EG92p . This finish the proof of item (3) of Theorem O □ 
At this moment, the proofs of item (4) and (5) are similar to the proofs of Corol- 
lary 4.3 and Theorem 4.4 in |DP05| respectively. The necessary minor modifications 
are similar to the ones treated in the proof of Lemma l3~4l and therefore we will skip 
the details. □ 
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In the same spirit of |AC81j . D. Daniclli and A. Petrosyan provided in |DP05 
a representation theorem for which Theorem 13.11 puts our minimizcrs u £ under the 
hypothesis of. The next Theorem will play an important hole in the investigation 
of fine regularity properties of the free boundary. 



Theorem 3.5. Let u e be a minimizer of problem \2,Hfy . Then 

(1) W" -1 (/C n d{u e > 0}) < oo for every compact set K C D c . 

(2) There exists a Borel function q s such that 

A pUe = fe H n - 1 L9K>o}, 

that is, for any ip G C^°(D ), there holds 

\Vu e \ p ~ 2 X7u E ■ X7ipdx = I i>q e dH n - Y . 

DC J{u c >0} 

(3) For any compact set K, C D , there exist positive constants c, C such that 

c<q £ <C 

and 

cr n-i < n n - l (B r {x)nd{u s > 0}) < Cr n -\ 



for every ball B x (r) C /C with x G d{u E > 0}. 
(4) For Ti."^ 1 almost all points in d{u e > 0}, an outward normal v = v(x) i 



defined and furthermore 

Ue(x +y) = q s {x)(y ■ v) + + o(y), 
where —* as \y\ — * 0. This allows us to define q e {x) = (u e ) u (x) at 



those points. 
(5) H n ~\d{u e > 0} \ d ied {u e > 0}) = 0. 

4. A GEOMETRIC-MEASURE HADAMARD'S VARIATIONAL FORMULA AND THE 

FREE BOUNDARY CONDITION 

In this section we suggest a geometric-measure version of the well known 
Hadamard's variational formula (see |G98j ) to deduce the free boundary condition 
of Problem 12. 21 Roughly speaking, given two points in the reduced free boundary, 
say Xi and X2, the idea is to make an inward perturbation around x\, an outward 
perturbation around Xi in such a way that we do not disturb very much the original 
volume and then compare the optimal configuration to the perturbed one in terms 
of the functional J e . Here are the details. 

Let p: K — » M be a nonnegative C°° function supported in [0,1], with, say 
J p(t)dt = 1. Hereafter, we shall write u = u e and fix two points x\ and X2 in the 
reduced free boundary d TC d{u > 0}. For any < r < dlst ^ X2 '> ; an d A > 0, we 
consider the vector field 



(4.1) P r {x) := ^ 



x + Xrp I j v{x\) x G B r (xi) 

x - Xrp ( — —J v(x 2 ) x G B r (x 2 ) 



r 

x elsewhere. 
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If v is any vector in R n , from direct computation, we see that 
(4.2) DP r (x) -v = v + (-iy +1 \\p' ( t^A) v ^ Xi) in Br{xi) 



\x - Xi\ 

Notice that, if A is small enough, P r is a diffeomorphism that maps B r {xi) onto 
itself. Indeed, if Asupp'(t) < 1, P r is a local injective diffeomorphism. Now, if 

[0,1] 

Ajo(t) < 1 - 1, for < t < 1, 

\P r {x) -Xi\<\x-Xi\ + Xrp (1^1) < r, 



= Id, therefore P r has to be 



for any x £ B r (xi). Finally, notice that P, 
onto. 

For each r > small enough, we will consider the r-perturbed configuration, v r 
implicitly defined by 

(4.3) v r (P r (x)) = u(x). 

The idea is to compare our optimal configuration {u > 0} to its perturbation 
{v r > 0} in terms of the penalized problem 12.21 An important geometric measure 
information we shall use is the blow-up limit. For any r > small enough and 
i = 1,2, consider the blow-up sequence, u l r : i?i(0) — > R, given by 

K(v) -u(xi + ry). 

From the blow-up analysis, we know, the set B x H {w* > 0} approaches {y G 
Bi | {y,v(xi)) < 0}, as r — > 0. Let us compute the change on the volume of the 
perturbation. More specifically, make use of the Change of Variables Theorem, we 
obtain 



\{B r (xi) n {v r > 0} 



(4.4) 



r " ifl r (ii)n{i r >o} 
dy 

B 1 n{v r (x i +ry)>Q} 



/ det (DP r (xi + ry)) dy 

JSiH{m^>0} 



l + (-iy+^\p'(\y\)/yy( Xi )\dy, 

B 1 n{( 3 /,^(x i )><0} 

as r — > 0. Notice there exists a constant C(p) so that, for any unit vector v e §" _1 
there holds 

(4-5) C(p)= J P \\y\)(±,v}dy. 

Bin{(y, I v)<0} 

A similar computation shows that 

|ffl r (^)n{u>o})| , 

(4-6) ^ \ dy, 



B l n{(y,u{x i ))<0} 
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as r — > 0. Combining 1)4.4(1 . 1)4.5(1 and 1(4. 6(1 . we conclude 

(|{«r>0}|)-(|{«>0}|) 
(4.7) - L _ o, 

as r — > 0. From the Lipschitz continuity of the penalization f £ , we obtain 

(4-8) / £ ((|K > 0}|)) - / B (|({u > 0}|)) < - e o{r n ). 

Now we shall turn our attention to the differential of the perturbation on the p- 
Dirichlet integral. Initially we observe that 

1 



\Vu{x)\*dx = / \Vu\.{y)\*dy 
(4.9) ' Jn 

\y<(y)\ p dy, 

Sin{M^>o} 

once d lc d{u l r > 0} is smooth. Now, applying twice the Change of Variables Theo- 
rem, taking into account that P r maps B r {xi) diffeomorphically onto itself, 
(4.10) 

' ' |Viv(a;)| p da; = \ [ ^P^p-^x))- 1 ■ Vu^O*))!*^ 



B r (xi) Br{Xi) 

1 

j.71 



\DP r {y)- 1 ■ Vu(y)\P\det (DP r (y))\dy 

B, ( , 

\DP r (xi +TZ)- 1 ■ Vui(z)\ p \det (DP r { Xl +rz))\dz. 
B 1 n{«* >o} 

Now, from 14.2(1 . using the fact that for any matrix A, with |^4| < 1, we have 

(Id + Ay 1 =Id+ £(-l)M\ we have 
i=i 

(4.11) DPrixi + rzy^VuUz) = V<(z)-(l) i+1 A^[^(z,V<(z))^(x l ) + o(A). 

\ z \ 

On the other hand, 

(4.12) | dct (DP r ( Xi + rzj) | = 1 + (-l) l+1 A^M(z, V { Xi )). 

\ z \ 

Combining ((4.9(1 . 1(4.10(1 . ((4.11(1 and 1(4.12(1 . we obtain 
(4.13) 



— J \Vv r (x)\ p - \Vu{x)\ p dx = (-ir +i A J \Vu l r {z)\ p ^-^{z,v{xi))dz 

B r (xi) Bin{u«>0} 

+ (-l) l A / p|V<(z)| p - 2 ^M(z,V<(z))(V<(z),K^))^ + o(A). 
Sin{u» >o} 

Again, from the blow-up analysis (see |DP05| ). for each 5 > 0, we know 

V< -» q(xi)v(xi), 
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uniformly in B\ n {(y,u(xi)) < —8}. Therefore, by r-uniform Lipschitz continuity 
of u\ , we have 

(4-14) V< -> -q{x i )v{x i )xB l n{(yM(x % ))<o}, 

in LP(Bi). Letting r ->• in lETTS]) . we find 
1 



|Vu r (cc)| p - |Vu(a;)| p da; — ► 

(4.15) r n '(\7\) 

(-iy +1 (p-l)X(q( Xi )) P / P -^{z,v{x % ))dz + o{\). 



Notice that 



Bin{«j.>o} 



div(p(|z|)) = 9 ^\ z,v(xi)), 



Thus, from Divergence Theorem and the blow-up analysis, 

(4.16) J P -^( ZiV ( Xl ))dz^ J p(\z\)dH n - 1 (z) = c(p). 

Sin{u*>o} B 1 n{{z,v(x i ))=o} 

Putting l|4.15|l and (|4.16|1 together, we obtain 



(4.17) / \Vv r {x)\ p - \Vu(x)\ p dx = r n X(p - l)c(p)(q(x 1 ) p - q(x 2 ) p ) + r n o(X). 
Jd c 

From the minimality property of u, l|4.8|l and 14.17(1 . 

(4.18) < J e (v r ) - J e (u) < r n X(p - l)c(p)(q( Xl f - q(x 2 f) + r n o(X) + - £ o{r n ) 
Dividing H4.18fl by r™ and letting r — > we obtain 

(4.19) < X(p - l)c(p)(q( Xi r - q(x 2 ) P ) + 0(A) 

Now dividing (|4.19() by A, letting A — ► 0, and afterwards reversing the places of x\ 
and x 2l we finally obtain 

(4.20) q{x 1 )=q{x 2 ) 

Since x\ and x 2 were taking arbitrarily in d re d{u > 0}, we have proven 
Theorem 4.1. There exists a positive constant X £ such that 

q e = X e , Vie d red {u £ > 0}. 

It now follows from DP05 that for each e > fixed, the reduced free boundary 
died{u E > 0} is a C l a smooth surface. Real analyticity of the reduced free boundary 
is then a consequence of KNS78 . It is worthwhile to point out that, for n = 2, a 
small variant of the main result in |DP| assures full regularity of the free boundary 
{u E > 0}, as long as p > 2 — er, for some universal constant a. 

In |C87| . |C88| and |C89 , L. Caffarelli introduced and developed the, by now, 
well known notion of viscosity solution of a given free boundary problem (for the 
Laplacian operator) . Our final goal of this section is to establish the free boundary 
condition obtained in Theorem 14. II in the viscosity sense. This will be used as a 
geometric tool in the remaining sections. 
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Theorem 4.2 (Free boundary condition in the viscosity sense). Let xq £ d{u £ > 0} 
be a free boundary point. Suppose there exists a touching ball B, i.e., dB n d{u £ > 
0} = {a;o}; such that either B C {u £ > 0} or B C {u £ = 0}. Then 

(4.21) u e (x) = \ £ {x - x a , v) + + o(\x - x \), 

where X £ is the positive constant provided in Theorem \4-l\ and v is the unit normal 
vector to dB, pointing inward to {u £ > 0}. 

Proof. Let B = B r (£). We shall first deal with the hypothesis that B C {u £ > 0}. 
For notation convenience, we will omit the center of the ball B. With no loss of 
generality, we can assume xo — and v = e n . Let be the auxiliary function 
solving 

( div(|Vii £ |P- 2 L>e) = inA = B r \S^ 

(4.22) I 6 = 1 on 8B r/2 
[ 6 = on dB r 

By the nondegeneracy of the gradient Vu e , see jPP05j . 9 £ C 1,a (A) and since 
it vanishes on dB, we have for some constant C = C(n,p, e, r) 

Q(x) = Cx n + o(\x\). 

Consider B to be a ball centered at such that B CC B. We define 

6q = sup |m : u £ (x) > m6(i) in B n s| . 

and any for k > 1, 

0k = sup |m : u e (x) > m6(a;) in 2 _fe B n b\ . 

The sequence {^fc} fe>1 is increasing and bounded by a constant i<T e , since u £ is 
Lipschitz continuous. Let — sup fc 0fe. If we set 

:= C£ 

we have 

u £ {x) > 0Q(x) + o(\x\) = 0x N + o(\x\). 

We claim, 

(4.23) u e (x) =9x N + o(\x\). 

Indeed, assume there exists a sequence {xk)k>i such that \xu\ = rf. — f and for 
some <5o > 0, 

u £ (x k ) - 0Q(x k ) > S Q \x k \. 
From the definition of 0k, we can choose a subsequence kj, where 2~ 6kj < r^ . < 
2- 4fe ^ and 

(4.24) u e (x k:j ) - § kj e(x kj ) > S Q \x kj \. 
Now, the function 

a(x) = u(x) - 0k j Q(x) 

satisfies div (\Vu £ \ p ^ 2 Da) = and a > in 2~ k iB n B. Since u e is Lipschitz 
continuous, from (|4.24l) . we conclude 

u £ {x) - kj h(x) > CS \x kj | 
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in a fixed portion of dB rk . From the Poisson representation formula, we have, in, 
\ B r k] n B 

(4.25) u s (x) - 9 kj @(x) >CS \x k] \. 
Now for j large enough, we have 

u e (x) > (9 kj + S^)@(x) 

which is a contradiction with the definition of 6 k . . We have proven so far, u E (x) — 
Ox 7v + o(|ac|). It remains to show, 9 = X E . Consider then the blow-up sequence 

(4.26) u*(x) = —u E (x + p k x), 

Pk 

as pk — > 0. By Lipschitz continuity and l|4.23[l . we know converges uniform in 
compact subsets to 

uf (a?) = 9{x,v) + . 

However, by a small modification of a classical argument, see |AC81| or |Fr82| . we 
can show the limit of a blow-up sequence of u E is an absolute minimizer of J E in 
any ball. Since the free boundary of it^° is smooth, we obtain from Theorem 14. II 
9 = A E , as desired. 

The case B r C {u e = 0}, follows by a small modification of the above arguments 
with the fact that Lemma Al in C89 also holds for nonnegative Lipschitz functions 
v satisfying div (\Vu e \p- 2 Dv) < 0. □ 

5. Recovering the original problem 

In this section we shall relate a solution to the penalized problem to a (possible) 
solution to our original problem. Roughly speaking the idea is that the function 
f E will charge a lot for those configurations that have a volume bigger than 1. Wc 
hope if the charge is too big, i.e., if e > is small enough, optimal configurations of 
Problem 12 . 21 will rather prefer to have volume 1 than paying for the penalization. 
We will follow the lines of |Teix05| . 

Lemma 5.1. There exist positive constants c and C , independent of e, such that 

c < |K > 0}| < 1 + Ce 

Proof. Let D* be any smooth domain containing D, so that \D* \ D\ = 1. From 
the minimality of u e , we have 

(5.1) J e (0= / \Vu e (x)\ p + f E (\{u E >0}\)<J E (u*)=C, 

J DC 

where u* is the p-harmonic function in D* \ D taking boundary data equal to tp on 
dD and on dD*. Thus 

i(|K>0}|-l) </ e (|K>0}|) <C. 

This proves the estimate from above. Let us turn our attention to the estimate by 
below. Expression l|5.1|l . together with Poincare inequality, provides 

(5.2) / \X7u E {x)\ p + \u E (x)\ p dx < C, 

J D c 

for some C independent of e. Let D$ be a tubular neighborhood of dD. For each 
x e dD fixed, let us consider the fiber F Xo :— {xo + tp(xo) | < t < 8} and denote 
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Q Xq :— Ti}{{u £ > 0} n F Xo ). From Mean Value Inequality, followed by Holder 
Inequality and l|5.2|) we have 



(5.3) 



Sip(xo) < / u(xa + tfi(xo))dt + / |Vtt(a;o + ifi(xo)) \tdt 



< 



cqILhv* i + 



s 



g + i 

Now we integrate l|5.3[) over 3D and obtain 

(5.4) / tpdS <C(S)\{u E >0}nD s \ 1/q . 

JdD 

Finally, from i|5.4[l . there must exist a constant, independent of e, so that \{u £ > 
0}| > c, as claimed. □ 

Lemma 5.2. There exists a positive constant C independent of e, so that X £ < C, 
where X £ is the constant provided by Theorem \4-l\ 

Proof. Applying Divergence Theorem to the field F\ = m|Vm| p_2 Vm, we have 

(5.5) / \Vu £ \ p dx= [ ip\Vu\ p - 2 d„udS, 

JD C JdD 

where [i is the outward unit vector in dD. If we apply Divergence Theorem to the 
field F 2 = |Vu| p ~ 2 Vu, we obtain 

(5.6) \ p e - 1 n n - 1 (d{u e >0})= [ iVuf^d^udS. 

JdD 

Isoperimetric inequality gives a universal bound by below to Ti n ~ 1 (d{u £ > 0}), i.e, 
H n ~ 1 (d{u £ > 0}) > c, for some c independent of e. Combining this with 15. 5|) and 
(|5.6() . we obtain 

\ £ <C(D,cp), 

as claimed. □ 

Lemma 5.3. There exists a universal positive constant c > 0, such that X £ > c, 
for all e > 0. 

Proof. Let zi £ D c be such that u £ (z 1 ) > for all e > 0. Let us denote by S = 
dist(,Zi, dD). Consider the smooth family of domains T t :— S j s +t (z 1 ) p| D c . Let t £ 
denote the first t such that T t touches d{u £ > 0}. Let us call xo = dT tc f]d{u £ > 
0}. Define ty £ to be p-harmonic function in T te \ To, with the following boundary 
values data: 

= mines and 9 S L~ =0. 

By the maximum principle we have u £ > \& £ in T te \ To- From Hopf's Lemma (see 
for instance |V84p we also know there exists a constant c > depending on 3D 
and inf ip, but independent of e, such that 

(5.7) *-v(*o)>c, 

where v denotes the outward unit normal vector of B± +t ^(zi), at xq. Recall from 
Theorem 14.21 we have the following asymptotic development around xq 

(5.8) ^(x) < u(x) =\ £ (x- x , v) + + o(\x - x \). 
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Dividing Ij5.8(l by \x — xo\, letting x — > xq and taking into account (|5.7|) . we finally 
obtain 

c < A e , 

as desired. □ 

We are ready to show the main theorem of this section. 

Theorem 5.4. If ' e is small enough, then any solution of Problem is a solution 
of Problem fOj) . 

Proof. Let us initially suppose \{u £ > 0}| > 1. In the same spirit of Section 0] 
consider a inward perturbation of the set {u £ > 0} with volume change V , in such 
a way that the set of positivity of the new function, u £ is still bigger than 1. Thus 

(5-9) / e (|{« e > 0}|) - / E (|K > 0}|) = --V. 

e 

From (|4.13|) and Lemma T5. 21 we have 

(5.10) jf o lV*r-|V«.|' = W + o(V) 

< C p V + o(V). 

Using the fact that J e (u £ ) < J £ (u £ ), l|5.9l) and l|5.1U|l . we find 

(5.11) < C p V + o(V) - -V. 

e 

Finally if we divide inequality (|5.11|) by V and let V — > 0, we obtain 

1 

£> c-v 

If \{u £ > 0}| < 1, we argue similarly, making an outward perturbation and using 
Lemma 15.31 to obtain another lower bound for e. Thus, if e is small enough, \{u £ > 
0}| automatically adjusts to be equal to 1. □ 

6. Radial Symmetry 

In this section we show a simple symmetry result of Problem (|1.2fl . Indeed, 
we shall show the best way of insulating a uniformly heated spherical body is by 
a ball. Recall, when (p = Constant, Problem l|1.2|) is equivalent to our original 
physical optimization problem. Here is the theorem: 

Theorem 6.1. Let D be the unit ball and ip = 1. Then Problem has a unique 

solution and it is radially symmetric. In particular the free boundary is a sphere. 

Proof. Let u ~ u £ be a solution to Problem 11.2f> . with D = B\ and ip = 1. Denote 
fi = {u > 0}. Let B ri and B r2 be the biggest ball inside f2 \ D and the smallest 
ball outside f2, respectively. Let yi £ dB ri C\dtt and y 2 E dB T2 Hdfl. Consider, hi, 
i = 1, 2 solution of 

!A p hi = in B ri \ B x 
hi = 1 on dBx 
hi = on dB Ti . 

It is simple to show hi is radially symmetric. Indeed, hi is the unique minimizer of 



E p (f) := f \Vf(x)\Pdx, 

JB ri \B! 
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among all functions / € W l p satisfying the according boundary data. For any 
orthonormal transformation O £ 0(n), consider 

h?(x) := hi(Ox). 
Clearly, hf has the same boundary data as hi and furthermore, 

E P (hf)= ( \0 T \7h l (Ox)\ p dx = E p (h l ). 

JB n \B! 

Thus hf = hi. Since O 6 0(n) was taken arbitrarily, hi has to be radial. 

In particular, the inward normal derivative of hi over dB ri is a positive constant 
Aj. Since r\ < we have 

(6.2) Ax > A 2 . 
Now, from the maximum principle, 

h\ < u < ft, 2 . 

Hence, from the free boundary condition in the viscosity sense, Theorem 14.21 we 
obtain 

(6.3) A 2 = {h 2 ) v {y 2 ) <X E < (/ii^G/i) - Ai. 

Combining Ij6.2|) and (|6.3(l . we conclude Ai — A 2 , and therefore, r\ — r 2 . This 
implies dfl has to be a sphere of radius r\ = r 2 . 

We have proven any solution to Problem (|1.2J) . with D = B\ and (p = 1 is 
radially symmetric. Uniqueness now follows due to the volume constraint. □ 
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